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• y—[ Abstract: High-permittivity dielectric particles with resonant magnetic 

p> I properties are being explored as constitutive elements of new metamaterials 

O and devices in the microwave regime. Magnetic properties of low-loss 

C/3 dielectric nanoparticles in the visible or infrared are not expected due to 

^ ^ intrinsic low refractive index of optical materials in these regimes. Here we 

C/2 analyze the dipolar electric and magnetic response of loss-less dielectric 

spheres made of moderate permittivity materials. For low material refractive 
index (< 3) there are no sharp resonances due to strong overlapping between 
different multipole contributions. However, we find that Silicon particles 
with refractive index ^3.5 and radius ~ 200nm present a dipolar and strong 
magnetic resonant response in telecom and near-infrared frequencies, (i.e. 
\j^ at wavelengths « 1.2 — 2/im). Moreover, the light scattered by these Si 

■t:;j- particles can be perfectly described by dipolar electric and magnetic fields, 

■^ quadrupolar and higher order contributions being negligible. 

in 

• © 2010 Optical Society of America 

^^ OCIS codes: (290.5850) Scattering, particles ; (160.1190) Anisotropic optical materials ; 

f^ (160.3820) Magneto-optical materials. 

> 

• ^H References and links 

y\: 1. H. C. van de Hulst, Light Scattering by small particles (Dover, New York, 1981). 

$—1 2. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley &Sons. New 

Cd York, 1998). 

3. M. I. Mishchenko, L. D. Travis and A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles 
(Cambridge Univ. Press, 2002). 

4. L. Novotny and B. Hecht, Principles ofNano-Optics, (Cambridge University Press, Cambridge, 2006). 

5. E. M. Purcell and C. R. Pennypacker, "Scattering and Absorption of Light by Nonspherical Dielectric Grains," 
Astrophys. J. 186, 705-714 (1973). 

6. B. T. Draine, "The discrete dipole approximation and its application to interstellar graphite grains," Astrophys. 
J. 333, 848-872 (1988). 



7. W. L. Barnes, A. Dereux, and T. W. Ebbesen, "Surface plasmon subwavelength optics," Nature 424, 824-830 
(2003). 

8. N. Engheta and R.W. Ziolkowski, "A Positive Future for Double-Negative Metamaterials," IEEE Trans. Microw. 
Theory Tech. 53, 1535-1556 (2005). 

9. J.B. Pendry, "Beyond metamaterials,"Nature Mater. 5, 763 (2006). 

10. V. M. Shalaev, "Optical negative-index metamaterials," Nature Photon. 1, 41-48 (2007). 

11. P. J. Schuck, D. P. Fromm, A. Sundaramurthy, G. S. Kino, and W. E. Moerner, "Improving the mismatch between 
light and nanoscale objects with gold bowtie nanoantennas," Phys. Rev. Lett. 94, 017402 (2005). 

12. F. Neubrech, A. Pucci, T. W. Cornelius, S. Karim, A. Garca-Etxarri, and J. Aizpurua, "Resonant plasmonic and 
vibrational coupling in a tailored nanoantenna for infrared detection," Phys. Rev. Lett. 101 , 157403 (2008). 

13. T. H. Taminiau, F. D. Stefani, F. B. Segerink, and N. F. Van Hulst, "Optical antennas direct single-molecule 
emission," Nat. Photonics 2, 234237 (2008). 

14. A. Alii and N. Engheta, "The quest for magnetic plasmons at optical frequencies," Opt. Express 17, 5723-5730 
(2009). 

15. A. Alii, and N. Engheta, "Dynamical theory of artificial optical magnetism produced by rings of plasmonic 
nanoparticles," Phys. Rev B 78, 0851 12 (2008). 

16. K. C. Huang, M. L. Povinelli, and J. D. Joannopoulos, "Negative effective permeability in polaritonic photonic 
crystals," Appl. Phys. Lett. 85, 543 (2004). 

17. M. S. Wheeler, J. S. Aitchison, J. I. L. Chen, G.A. Ozin, and M. Mojahedi, "Infrared magnetic response in a 
random silicon carbide micropowder," Phys. Rev. B 79, 073103 (2009). 

18. L. Jyhla, I. Kolmakov, S. Maslovski and S. Tretyakov, "Modeling of isotropic backward-wave materials com- 
posed of resonant spheres," J. Appl. Phys. 99, 043102 (2006). 

19. L. Peng, L. Ran, H. Chen, H. Zhang, J. A. Kong, and T. M. Grzegorczyk, "Experimental Observation of Left- 
Handed Behavior in an An'ay of Standard Dielectric Resonators," Phys. Rev. Lett. 98, 157403 (2007). 

20. J. A. Schuller et al., "Dielectric Metamaterials Based on Electric and Magnetic Resonances of Silicon Carbide 
Particles," Phys. Rev. Lett. 99, 107401 (2007). 

21. K. Vynck, D. Felbacq, E. Centeno, A. I. Cabuz, D. Cassagne, and B. Guizal, "All-Dielectric Rod-Type Metama- 
terials at Optical Frequencies," Phys. Rev Lett. 102, 133901 (2009). 

22. R. K. Mongia, and P. Bhartia, "Dielectric resonator antennas-A review and general design relations for resonant 
frequency and bandwidth," Int. J. Microwave Millimeter-Wave Comput. -Aided Eng. 4, 230247 (1994). 

23. R. C. J. Hsu, A. Ayazi, B. Houshmand, and B. Jalali, "All-dielectric photonic-assisted radio front-end technol- 
ogy," Nat. Photonics 1, 535538 (2007). 

24. Jon A. Schuller and Mark L. Brongersma, "General properties of dielectric optical antennas," Opt. Express 17, 
24084-24095 (2009). 

25. Q. Zhao, J. Zhoul, F. Zhang, D. Lippens, "Mie resonance-based dielectric metamaterials," Materials Today 12, 
60-69 (2009). 

26. E.D.Palik, Handbook of Optical Constants of Solids (Academic Press, Orlando, Florida, 1985). 

27. S. Albaladejo,R. Gomez-Medina, L. S. Froufe-Perez, H. Marinchio, R. Carminati, J. F. Toirado, G. Armelles, 
A. Garcia-Marti'n and J.J. Saenz, "Radiative corrections to the polarizability tensor of an electrically small 
anisotropic dielectric particle," Opt. Express 18, 3556 (2010). 

28. PD. Garcia, R. Sapienza, A. Blanco, and C. Lopez, "Photonic Glass: A Novel Random Material for Light," Adv. 
Mater. 19, 2597-2602 (2007). 

29. M. Reufer, L.F. Rojas-Ochoa, S. Eiden, J.J. Saenz and F. Scheffold, "Transport of Light in Amorphous Photonic 
Materials," Appl. Phys. Lett. 91, 171904 (2007). 

30. M. Ibisate, D. Golmayo and C. Lopez, "Photonic Ciystals: Silicon Direct Opals," Adv. Mater. 28, 2899 - 2902 
(2009) 

31. S. Albaladejo, M.I. Marques, M. Laroche and J.J. Saenz "Scattering Forces from the Curl of the Spin Angular 
Momentum of a Light Field," Phys. Rev. Lett. 102, 1 13602 (2009). 

32. PC. Chaumet and A. Rahmani, "Electromagnetic force and torque on magnetic and negative-index scatterers," 
Opt. Express 17, 2224-2234 (2009). 

33. M. Nieto-Vesperinas, J.J. Saenz, R. Gomez-Medina and L. Chantada, "Optical forces on small magnetodielectric 
magnetic particles," Opt. Express 18, 1 1428 (2010). 



1. Introduction 

Electromagnetic scattering from nanometer-scale objects has long been a topic of great interest 
and relevance to fields from astrophysics or meteorology to biophysics and material science 
Hill] [3] mis] IS- During the last decade nano-optics has developed itself as a very active field 
within the nanotechnology community. Much of it has to do with plasmon (propagating) based 
subwavelength optics and applications |7], the synthesis of negative-index optical metamateri- 







Fig. 1. Incident field vector 



als iSEIIOl or the design of optical antennas imi fT2l[T3l . 

Magnetic effects, a key ingredient of relevant microwave applications, cannot be easily ex- 
ploited in the optical range (visible or infrared) due to intrinsic natural limitations of optical 
materials. The quest for magnetic plasmons and magnetic resonant structures at optical fre- 
quencies 1 14 j has been mainly focused on metallic structures with the unavoidable problems of 
losses and saturation effects inherent to metals in the optical range. Moreover, their size is often 
required to be comparable to, or a fraction of, the operating wavelength in order to provide a 
non-negligible response, and their magnetic response is strongly affected and influenced by an 
intrinsic quadrupolar contribution, which introduces additional radiation losses and impurity of 
radiation and polarization 1 15 1. 

As an alternative approach, there is a growing interest in the theoretical and experimen- 
tal study of high-permittivity dielectric objects as constitutive elements of new metamaterials 
QH [nl [I8l [m 1^ EH and antennas based on dielectric resonators |!22l |23l |24| . Interesting 
magnetodielectric properties are usually associated to low loss and large dielectric constants 
which, being accessible at microwave and terahertz frequencies, are still a challenging issue in 
the infrared (IR) and visible frequency ranges |25|. However, it has recently been shown 1211 
that silicon rod arrays could exhibit a true metamaterial left handed dispersion branch in the 
visible to mid-IR despite the moderate value of the refractive index {m ^ 3.5). 

Motivated by these results, we have analyzed in detail the magnetodielectric properties of 
loss-less dielectric spheres made of moderate permittivity materials in the infrared regime. Al- 
though the analysis of scattering resonances of homogeneous spheres have been analyzed in 
the past HlllKSj, we concentrate here on the relevance of the first magnetic dipolar resonance 
in non-absorbing dielectric nano-spheres (e.g. silicon spheres in the mid-IR). As we will see, 
submicron Silicon spheres (with radius a ~ 200nm) can present a strong magnetic resonance 
response at telecom and infrared frequencies, (namely, wavelengths X w 2na w 1.5/im). More- 
over, based on standard Mie theory [!_], we have found that the scattering of these Si nanoparti- 
cles can be perfectly described by dipolar electric and magnetic fields; quadrupolar and higher 
order contributions being negligible. Interestingly, the magnetic resonance line remains well 
defined and resolved even for lower index materials, as Rutile-like Ti02 with an average index 
m w 2.8 in the near-IR, even though in this latter case the electric dipole resonance line strongly 
overlaps with the magnetic dipole one, as well as with those of higher order multipoles. 

2. Extinction resonances of a dielectric sphere. The Mie theory revisited 

Consider a non-absorbing dielectric sphere of radius a, index of refraction nip and dielectric 
permittivity Ep = m^ in an otherwise uniform medium with real relative permittivity £/, and re- 



fractive index nif, — y/eij. The magnetic permittivity of the sphere and the surrounding medium 
is assumed to be 1. Under plane wave illumination, and assuming linearly polarized light, the 
incident wave is described by 

E = Eouze'^e-'"" , B = BoUre^'^e-'"" (1) 

where k — mi,co/c = m/,2;r/A, A being the wavelength in vacuum and Bo = Mo^o = ^ {nVi/c)EQ 
(see Fig.fT]). The field scattered by the sphere can be decomposed into a multipole series, (the 
so-called Mie's expansion), characterized by the electric and magnetic Mie coefficients {«„} 
and {b„}, respectively; (ai and b\ being proportional to the electric and magnetic dipoles, a2 
and b2 to the quadrupoles, and so on). We shall find useful to write the Mie coefficients in terms 
of the scattering phase-shifts a„ and j3„ UJ 

a„ = -(l-e"2'"«) =/sina„e"'"" (2) 

bn - \(l-e-'-'^")^i^m^ne-'^-^ (3) 
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where 



j^jj^ ^ m^injy) [xin{x)\ - in{x) \yin{y)'\ ,4, 

ni^jn {y) [xyn (jc)] ' - yn {x) [yjn (y)] " 

^^n ^ jn {y) [xjn {xy\ - ;„ jx) [yjn (j)]' 

jn{y) [xy„{x)^ -yn{x) [yj„(yy\ ' 

m = nip/mij being the relative refractive index, x — ka the size parameter and y = mx. j„{x) 
and yn{x) stands for the spherical Bessel and Neumann functions, respectively, and the primes 
indicate differentiation with respect to the argument. In absence of absorption, i.e. for real m, 
the phase angles a„ and )3„ are real, then the extinction and scattering cross sections, ffext and 
as, respectively, have the common value (H] 

ffS = C^ext = T2" L (2n + 1) {sin^ «„ + Sin^ Pn} = Y, i^E,n + OM,n} (6) 

'^ n=l n=l 

In the small particle limit (x <C 1) and large particle permittivities (m 3> 1) the extinction 
cross section presents characteristic sharp resonance peaks. The values of y — mx at which 
the angles a„ or j3„ are n/2,3n/2, . . . ,etc, define the resonance points. At each resonance, the 
extinction cross section is of the order of A^ and it is independent of either the particle size or 
the refractive index yj 

^r=C^SI{.«,;,„»l}«^(2«+l) (7) 

Asymptotically, the first resonance peak occurs at y = TT (i.e. A = m2a) corresponding to the 
magnetic dipole term of coefficient bi. 

Interesting consequences and applications of well defined Mie resonance lines, associated 
to low loss and large dielectric constants, are accessible for different materials at microwave 
and terahertz frequencies. However, as m decreases there is an increasing overlap between the 
wavelength dependent cross-section peaks, and the sphere resonant response smears out. Since 
usually non-absorbing materials present low refractive index in the infrared (IR) and visible 
frequency ranges, Mie resonances of small particles in these regimes have not been considered 
in detail. However, the above analysis on the cross-section, carried out for Si nanoparticles 




Fig. 2. Scattering cross section map of a non-absorbing Mie sphere as a function of the 
refractive index m and the size parameters = 27ta/X. Green areas correspond to parameter 
ranges where the magnetic dipole contribution dominates the total scattering cross section, 
while red areas represent regions where the electric dipole contribution is dominating. The 
remaining blue-saturated areas are dominated by higher order multipoles. Brightness in the 
color-map is proportional to the total cross section. White horizontal lines represent the 
j-range covered by figures[3](m = 3.5) and|5](m = 2.8) 



in the IR, indicate that the aforementioned asymptotic behavior can be extended to a relative 
refractive index as low as m w 3. 

In Fig. |2] we show the scattering cross section map of a non-absorbing Mie sphere as a 
function of the refractive index m and the size parameter x = 2na/X. Brightness in the color- 
map is proportional to (Text while color code represents the contribution of electric and mag- 
netic dipoles to the total cross section. The RGB (Red, Green, Blue) code is formed by taking 

RGB — OeaR + <7m,iG + O^^sB with dres = y'^ixt ~ '^e i^ Ca/.i- The whole map is normal- 
ized to avoid over-saturation. Hence, green areas correspond to parameter ranges where the 
magnetic dipole contribution dominates the total scattering cross section, while red areas repre- 
sent regions where the electric dipole contribution is dominating. The remaining blue-saturated 
areas are dominated by higher order multipoles. 

In the micrometer wavelength regime, within the transparent region of silicon (A > 1 100 
nm), the index of refraction can well be approximated by a real constant nip « a/12 ^3.5 
(see for example 1261 ). The scattering (or extinction) cross section of a Si sphere of radius 
a = 23Qnm in vacuum (ot/, = 1) is plotted in Fig. [3] Although there is a partial overlap between 
the first dipolar peaks, the magnetic line of the first resonant peak (at A « 2ma) is still very well 
resolved. Interestingly, for wavelengths larger than A w 1200nm, the cross section is completely 
determined by the first bi and ai coefficients. In other words, in this regime Si particles can be 
treated as non-Rayleigh dipolar particles. 
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Fig. 3. Scattering cross-section a^ versus tiie wavelengtli A for a 230nm Si spliere (the 
refraction index m = 3.5 is constant and real in tliis wavelengtfi range). Tlie contribution of 
each term in the Mie expansion is also shown. The green line corresponds to the magnetic 
dipole contribution. 



3. Magnetic and electric resonances of dipolar particles 

Let us consider further in some detail the scattering properties of small dielectric particles 
(x <C 1). In the .jc <C 1 limit, the sphere is sufficiently small so that only the dipole scattered 
fields are excited. The induced dipole moments, proportional to to the external (polarizing) 
fields, E and B, are usually written in terms of the particle electric and magnetic polarizabilities 
Ue and Um, respectively : 



p = eoe/,o:£E 



m = £o£/,aMB; 



where 
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By using the definition of the phase angles cci and j3i, [cf. Eqs. (J2| - (J5]l], we can rewrite the 
polarizabilities as 
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Fig. 4. Maps for the square modulus of the total electric an magnetic vectors, emitted by 
the Si nanoparticle of radius a = 230nm under plane wave illumination, (cf. Fig. [TJ. The 
left and right panels correspond to the wavelengths I250nm and 1680nm of the electric 
and magnetic resonance peaks of Fig. [3] respectively. The corresponding far field radiation 
patterns are shown in the bottom. 
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Fig. 5. Scattering cross-section <Js versus the wavelength A for a sphere of radius a = 
150«m (its relative refractive index is m = \/8 = 2.82, constant and real in this wavelength 
range). The contribution of each coefficient in the Mie expansion is also shown. The green 
line corresponds to the magnetic dipole contribution. 



are the quasi-static polarizabilities. In terms of these magnitudes, the extinction and scattering 
cross sections then read 



ffs 



= klm{aE 



am} 



f) 



(12) 
(13) 



Of course, in absence of absorption, a^ and a]^ are real quantities and we recover the well 
known optical theorem result (Text = O5. 

In particular, in the Rayleigh limit, when y . 
electrostatic form. 



2nma/X ^1, «£ and a]^ approach the 
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From Eq. ( 14 1, together with Eq. (fTO|, we recover the well known expression for the polariz- 
ability of a Rayleigh particle including radiative corrections [6, 27J. Notice that in this limit 
the magnetic polarizability is negligible. Very small particles always behave as point electric 
dipoles. 

However, asy increases, (i.e. as X decreases), there is a crossover from electric to magnetic 
behavior as shown by Fig. plfor a Si particle. The y values at which the quasi-static polariz- 
ability ajj ( Ctj^ ) diverges define the electric (magnetic) dipolar resonances, f^ear the first 



bi-resonance, the particle essentially behaves like a magnetic dipole(cf. Fig. 2). If A decreases 
further, ai peak dominates and the sphere becomes again an electric dipole. Notice however 
that, due to the overlap between the electric and magnetic responses, the radiation field near the 
resonances does not correspond to a fully pure electric or a fully pure magnetic dipolar exci- 
tation. Yet, as seen in Fig. 2 for the aforementioned Si sphere, the magnetic dipole contribution 
is, at its peak, about five times larger than the electric dipole one. 

As a further example, we have plotted in Fig. |4] both the near field intensity maps and far 
field radiation patterns at the resonant peaks of ai and bi, (cf. Fig.l3]l, numerically calculated 
from the full Mie solution (and not only from the dipolar approximation). The slightly distorted 
patterns with respect to the single dipole case are completely explained in terms of the sum of 
these two mutually coherent electric and magnetic dipoles. 

Silicon nanoparticles with radius of the order of 200-250nm are then fully characterized by 
electric and magnetic polarizabilities as given by Eqs. [TOJin the mid and near-IR region. By 
tuning the appropriate resonance wavelength. Si particles exhibit a strong magnetic response, 
(only slightly perturbed by their electric dipole counterpart). 

It should be remarked that if the relative refractive index is smaller, the overlap between the 
Mie resonance lines increases and then, a particular one single Mie term magnetic response is 
not well defined. As an example, in Fig.|5]we represent the scattering cross section of a 150nm 
radius particle with m — \/8 « 2.8 (which would be of the order of the averaged refractive index 
of Ti02 (rutile) in the near IR ll26l ). As seen in this figure, although beyond 700nm the scatte- 
ring is still described by the sum of the electric and magnetic dipoles, the peak corresponding 
to the electric dipole, ai, now disappears, and the overlapping between lines that produces this, 
manifests a significant contribution of this electric dipole at the magnetic resonance wavelength. 

4. Conclusion 

We have analyzed the dipolar electric and magnetic response of lossless dielectric spheres made 
of both low and moderate permittivity materials. Based on the Mie expansion, we have derived 
general expressions for the electric and magnetic polarizabilities of dielectric spheres. We found 
that submicron Silicon particles present a strong magnetic resonant response in the mid-near 
infrared. Interestingly, the light scattered by these Si nanoparticles of appropriate size is per- 
fectly described by dipolar electric and magnetic fields, being quadrupolar and higher order 
contributions negligible. These results can play an important role, not only in the field of meta- 
materials or optical antennas, but also in tailoring the light transport through complex dielectric 
media like photonic glasses ll28l l29ll with intriguing magnetic properties. Fabrication of new 
materials made of highly monodisperse subwavelength silicon spheres 1301 may then lead to 
a new generation of magnetodielectric optical materials. At the magnetic or electric resonance 
wavelengths the extinction cross section is of the order of A^ reaching its maximum theoretical 
limit (independent of the particle size or refractive index) |[T] ( see also ll24l ). The large dipolar 
cross-section of magnetodielectric spheres near a magnetic resonance would also imply strong 
radiation pressure magnetic forces 113111321 [33l leading to new concepts related to the optical 
forces on magnetic particles . 
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